Scaling of the thermal resistivity of ^He in restricted geometries 
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The thermal resistivity and its scaling function in quasi-2D *He systems are studied by Monte 
Carlo and spin-dynamics simulation. We use the classical 3D XY model on L x L x H lattices 
with L ^ H, applying open boundary condition along the H direction and periodic boundary 
conditions along the L directions. A hybrid Monte Carlo algorithm is adopted to efficiently deal 
with the critical slowing down and to produce initial states for time integration. Fourth-order 
Suzuki- Trotter decomposition method of exponential operators is used to solve numerically the 
coupled equations of motion for each spin. The thermal conductivity is calculated by a dynamic 
current-current correlation function. Our results show the validity of the finite-size scaling theory, 
and the calculated scaling function agrees well with the available experimental results for slabs using 
the temperature scale and the thermal resistivity scale as free fitting parameters. 

PACS numbers: 64.60Ht, 67.40Kh, OS.lOLn 
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I. INTRODUCTION 

High-resolution measurements of various prop- 
erties of ^He confined in restricted geometries 
near the bulk superfluid transition temperature Tx 
have been extensively carried out for over three 
decadeai.^.3,4,5,6,7,8,9,io,ii,i2,i3,i4_ The measurements 

approach so close to T\ that the correlation length 
becomes macroscopic in size. As a result, the whole 
fluid acts in a correlated way and the values of global 
properties are changed. This offers an opportunity for 
testing the finite-size scaling theory which describes the 
effect of confinement in a finite geometry near a critical 
point. 

Physical systems which exhibit a second-order phase 
transition and which are confined in a finite geometry 
(e.g., a film, a pore, or a box) are thought to be well 
described by the finite-size scaling theory^^ at tempera- 
tures close to the critical temperature T\ . The finite-size 
scaling theory is based on the idea that the finite-size 
effects can be observed when the correlation length ^ be- 
comes of the order of the finite system size (i.e., the side 
of the cube, the thickness of the film, or the diameter of 
the pore). For a physical quantity O, this statement can 
be expressed as followai^: 



0{t,H) 



0{t,H = oo) 



(1) 



where 



t = 



Tx 



and X = 



H 



Here H denotes the relevant confining length and 
^{t, H = oo) is the correlation length of the bulk sys- 
tem, /(x) is a universal function which does not depend 
on the microscopic details of the system. It does, how- 
ever, depend on the observable O, the geometry of the 
system, and boundary conditions imposed on the system. 



Liquid helium ''He has been widely used to examine 
the validity of the finite-size scaling theory of critical phe- 
nomena; sophisticated experimental studies were carried 
out in microgravity environment, for example, Lipa et 
ali^ measured the specific heat of helium confined in 
a parallel plate geometry with a spectacular nanokelvin 
resolution, and the shape of the confining geometry, such 
as film, pore, and box, can be designed with such a pre- 
cision that the relevant confining length is well defined. 

For static properties, the scaling of superfluid ^He has 
been studied experimentally, analytically and numeri- 
cally. For example, the speciflc heat near the superfluid 
transition of '^He has been measured for confinements 
which vary by a factor of over 1000 and the data of a large 
extent can be collapsed upon a unique function when 
properly reducediiSiii. Field-theoretical calculations for 
the standard Landau-Ginzbrug free energy functional in 
different geometries with dirichlet boundary conditions 
have been carried ouli^iSii and the results agree with the 
specific heat measurements^iS. Schultka and Manousakis 
calculated the superfluid density and the specific heat of 
^He in film and pore geometries using Monte-Carlo sim- 
ulation and they demonstrated the validity of the finite- 
size scaling theorjsi2iii*iS. 

Beside static properties, finite-size scaling theory can 
be also understood by studying dynamical and transport 
properties near the critical point. Among them, the ther- 
mal conductivity A of '*He is a good candidate, because 
it is a measurable property and its bulk transition has 
been carefully studiediSiiSiSSiSi. D. Murphy et al^ and 
E. Genio et alm^ measured the thermal conductivity of 
liquid ^He confined in a " microchannel plates" of thick- 
ness 2 mm with holes 0.5 /xm and 1 /im in diameter, and 
their data are consistent with a universal scaling function. 
Nho and Manousakisi^i studied the thermal conductivity 
A of ''He confined in a pore-like geometry and their simu- 
lational results agree well with the experimental results. 
Rather recently, Jerebets et al. measured the thermal 
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conductivity of liquid ^He confined in a film geometry24'. 
In this paper we will study the thermal conductivity A 
of confined helium using Monte Carlo and spin-dynamics 
methods, calculate the scaling function associated with A 
for a film geometry, and compare it with the experimen- 
tal results. 

The remainder of the paper is organized as follows: In 
Sec. II, we will present the model, the simulation meth- 
ods, and the method for extracting the thermal conduc- 
tivity. In Sec. Ill, we will examine the finite-size scaling 
theory for the thermal resistivity and compare our cal- 
culated scaling function with experimental results. We 
summarize our results in the last section. 



II. MODEL AND SIMULATION METHOD 

In this section, we will briefly describe the model and 
simulation methods used to study the properties of su- 
perfluid ^He and show how the thermal conductivity is 
computed in our model. Matsubara and MatsudaSSi have 
proposed a lattice model to explain the properties of liq- 
uid helium. In the model, the liquid is regarded as a 
lattice composed of atoms and holes. In terms of oper- 
ators which create or annihilate an atom at each lattice 
point, it is proved that the lattice model is equivalent to 
the spin system. Both systems are classified as model F 
(or E in the absence of an external field) in the classifi- 
cation of dynamical models and same universality class 
with the XY model, so we can use the XY model to study 
the properties of liquid ^ReS^iSL. 

The Hamiltonian of the XY model on a lattice takes 
the following form: 



ciated with the lattice point i is defined by 



-J ^ (SfS^ 
<ij> 



(2) 



where J sets the energy scale, < ij > denotes a nearest 
neighbor pair of spins on a simple lattice in three dimen- 
sions, and the summation is over all nearest neighbors. 
In this model, each spin is a classical spin with three 
components, St — {Sf , Sf , Sf). 

In our calculations, we use a film geometry, i.e. a 
L X L X H lattice with L ^ H, in order to mimic experi- 
ment. Open boundary conditions are imposed along the 
confining dimension (the H direction). In our model the 
spins at the open boundary have no neighbors outside the 
confining space. Periodic boundary conditions are used 
along the large planar dimensions ( L direction) , because 
they approximate the limit L oo better. 

The thermal conductivity A of liquid ''He at a given 
temperature T can be calculated using the dynamic 
current-current correlation functionSSiSS 



A oc 



(3) 



= Jisfst 



QX qy \ 



(4) 



where the notation i + denotes the nearest neighbor of 
the lattice site i in the z lattice direction. 

The thermal conductivity calculated by above equa- 
tions strongly depends on the initial configuration of the 
system. The thermal conductivity at a given temperature 
T is an average value for all possible initial configurations. 
Therefore, a sequence of uncorrelated equilibrium config- 
urations is needed to provide starting points for the spin 
dynamics. 

These configurations are obtained from a Monte Carlo 
simulation using the Hamiltonian given by Eq.(2). We 
use a hybrid Monte Carlo procedure^^, which consists 
of 10 updates (MCCMOCMCCO), where M is Metropo- 
lis updateSS, C is cluster update^ and O is the over- 
relaxation algorithm's. 

Using this hybrid algorithm, we generate approxi- 
mately 2000 uncorrelated configurations from the equi- 
librium canonical ensemble at a given temperature. From 
each configuration we do dynamic simulation according 
to the equations of motion for the XY model, which are 
given as follow8SSi2&: 



At as. 



(5) 



where the z component of the current density ji asso- 



Starting from a particular initial spin configuration, we 
perform numerical integration of these equations of mo- 
tion using a recently developed decomposition method^i, 
which guarantees exact energy conservation and conser- 
vation of spin length \Sk\ = 1 and conserves Mz within 
its numerical truncation errors. The integration is car- 
ried out to a maximum time tmax (typically of the order 
of tmax = 200) to make sure that we determine the real- 
time history of every configuration within a sufficiently 
long interval of time {0 < t < tmax)- The time step 
is Ai — 0.1, which guarantees sufficient accuracy with 
respect to the conservation of ■ Fourth-order integra- 
tors are used, which are more accurate than a second- 
order method for same time step. Finally, we compute 
the average of a time-dependent observable (such as the 
current-current correlation function) over all results rel- 
ative to all the independent initial equilibrium configu- 
rations generated via the hybrid Monte Carlo procedure. 
All error bars for the thermal resistivity correspond to 
one standard deviation. 



III. RESULTS 

In this section, we calculate the thermal resistivity, ex- 
amine its scaling behavior with respect to 77, and com- 
pare the scaling function with the experimental results. 
For a film geometry, L x L x iJ, it is better to keep 
H finite and let L — > oo such that the system can not 
feel the effect of L for any temperature. In other words. 
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FIG. 1: The thermal resistivity with a fixed = 6 at as 
a function of L/ H . 



we wish the correlation length to be larger than H , but 
always less than L for any temperature. We can realize 
this by performing extrapolation to the infinite limit from 
the results on finite L. However, we do not have to do 
this, because we find that by applying periodic bound- 
ary conditions along the directions of L, the finite-size 
effects due to L are already not significant as long as 
L > 5H (see Fig. 1) . Fig. I shows the thermal resisitiv- 
ity R{t,H) = l/X{t,H) for film lattices as a function of 
L/H at the critical temperature, where the correlation 
length is infinite. From Fig. I, we can see that within 
the error bars, the results of thermal resistivity are same 
for L/H > 5. 

Fig. 121 shows some of our results of the thermal re- 
sistivity -R(r, H) as a function of temperature T for 
various lattice sizes H. The bulk transition temper- 
ature T\ — 1.5518(2) is obtained from Monte Carlo 
simulation^^. At high temperature, there is no significant 
scaling behavior and the thermal resistivity for different 
H has no obvious difference, because at high tempera- 
tures, £,{t) < H. Near T\, the data show strong effects of 
confinement. The smallest size shows the greatest round- 
ing of transition and the thermal resistivity is biggest for 
smallest size at the critical temperature. 

Now we would like to check the finite-size scaling hy- 
pothesis for the thermal resistivity and compare our re- 
sults with the existing experimental results'^"'. The de- 
pendence upon t of the bulk thermal resistivity can be 
described by a power law 



R{t) = Rot^, 



(6) 



where x is the dynamic critical exponent of the ther- 
mal resistivity. The dependence upon t of the correlation 
length can be described by a power law 



(7) 



where u is the critical exponent of the correlation length. 
Using Eq.(l), the finite-size scaling expression for the 




FIG. 2: Thermal resistivity R{T, H) versus temperatures for 
film geometry with sizes that correspond to H =4, 6, 10 and 
L — 5H. The bulk critical temperature T\ — 1.5518 is also 
shown. 



thermal resistivity R{t, H) is given by 



R{t,H)H 



g{tH 



(8) 



where g{x) is a universal function"^'*. 

Fig. 3 shows a scaling plot of the thermal resistivity 
scaling function g{x) = R{t, H)H^/^ versus the scaled 
reduced temperature parameter x = tH^^^ , where the 
reduced temperature is taken relative to the bulk tran- 
sition temperature T\. We use the value v = 0.6705 as 
determined by Goldner and Ahlers^^^ and x — 0.44 deter- 
mined by Ahlers^^. From Fig. 3, we can find that near 
T\ our simulation data collapse onto one single curve for 
a wide range of values of H and i, thus supporting the 
validity of finite-size scaling theory. When T is large, 
where ^(t) < /f, we can not obtain the unviersal func- 
tion Eq.(8) using Eq.(l), therefore, the simulation data 
might not collapse onto one single curve. 

In Fig. 3 we also compare our universal function g{x) 
with the experimental data obtained by Jerebets et ali^ 
represented by a solid line. In order to do this, we used 
two multiplicative constants as free fitting parameters, 
one multiplying the scale of a;-axis and another the scale 
of y-axis. The agreement between simulation and exper- 
iment is quite good. 

The geometry is one of important factors which, in 
principle, can determine the scaling function. In Fig. 
4, we compare the scaling functions for films and pores. 
From Fig. 4 we can find that near T\ the scaling function 
for films and pores are different alghough their shapes 
are similar. Compared with scaling function for films, 
the scaling function for pores shifts to the left, further 
away from the bulk behavior. This is expected; since in 
comparison with the film geometry, the pore geometry 
restricts the dimensionality of the system more, thereby 
further limits critical fiuctuations. As a consequence the 
thermal resistivity for pore geometry is higher at Tx, and 
decreases to at lower temperature. 
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o 4x20x20 
□ 6x30x30 
A 10x50x50 
— S. Jerebets et al. 




IV. SUMMARY 



FIG. 3: The universal function g{x) obtained for film ge- 
ometry. The solid line represents the available experimen- 
tal results for film geometry. In the experimental results the 
resistivity scale and the temperature scale are used as free 
parameters. 




FIG. 4: Comparison of the finite-size scaling functions ^(a;) 
for films and pores. 



We have calculated the thermal resistivity R{t, H) of 
liquid ''He in a film geometry (on a L x L x H lattice) 
applying periodic boundary conditions along L directions 
and open boundary conditions in the H direction. We 
were able to demonstrate the validity of finite-size scal- 
ing theory and we obtained the thermal resistivity scal- 
ing function g{x) using known values for the critical ex- 
ponents and no adjustable parameters. We find good 
agreement for scaling functions between simulational and 
experimental results using the temperature scale and the 
thermal resistivity scale as free parameters. We also com- 
pared our calculated scaling function for film geometry 
with the result for pore geometry and found a systematic 
shift. 
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